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We present the analytical solutions for the evolution of matter density perturbations, for a model 
with a constant dark energy equation of state w but when the effects of the dark energy perturbations 
are properly taken into account. We consider two cases, the first when the sound speed of the 
perturbations is zero = 0 and the general case 0 < < 1. In the first case our solution is 

exact, while in the second case we found an approximate solution which works to better than 0.3% 
accuracy for k > WHo or equivalently k/h > 0.0033Mpc“^. We also estimate the corrections to the 
growth index 7 ( 2 ), commonly used to parametrize the growth-rate. We find that these corrections 
due to the DE perturbations affect the growth index 7 at the 3% level. We also compare our new 
expressions for the growth index with other expressions already present in the literature and we find 
that the latter are less accurate than the ones we propose here. Therefore, our analytical calculations 
are necessary as the theoretical predictions for the fundamental parameters to be constrained by the 
upcoming surveys need to be as accurate as possible, especially since we are entering in the precise 
cosmology era where parameters will be measured to the percent level. 

PACS numbers: 95.36.+X, 98.80.-k, 98.80.Es 


I. INTRODUCTION 


Under the assumptions that the Universe at large 
scales is homogeneous and isotropic, that it can be de¬ 
scribed by General Relativity or some other modified 
gravity theory, such as f{R) or any other metric theory 
whose effects at the perturbations level can be taken into 
account by the effective Newtonian constant Geff(a), see 
Refs. [IHl], and finally under the subhorizon approxima¬ 
tion {k 3> aH), then it can be shown that the growth of 
matter is governed by the second order differential equa¬ 
tion: 




3 Umj,Geff(a)/GAr^(a) 




= 0 , 
( 1 ) 
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where 6 = is the matter density contrast that de¬ 
scribes the growth of matter (known as the growth fac¬ 
tor), H{a) is the Hubble parameter, flmo is the matter 
density today and Hq is the Hubble constant. 

Making the assumption that the dark energy com¬ 
ponent can be described by a constant equation of 
state w and negligible dark energy perturbations, i.e. 
Ges{a)/GN = 1, then Eq. 0 can be easily solved an¬ 
alytically. The differential equation 0 has in general 
two solutions that correspond to two different physical 
modes, a decaying and a growing one, that in a mat¬ 
ter dominated Universe in GR behave as (5 = and 

as (5 = a respectively. Since we are only interested in 
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the latter, we demand that at early times a™ 1, usu¬ 
ally during matter domination, the initial conditions have 
to be chosen as S(ain) — am and S'(am) — 1. When 
Ges{a)/Gtq = 1 we get GR as a subcase, while in gen¬ 
eral for modified gravity theories, the term Geff can be 
scale and time dependent. 

For a flat GR model with a constant dark energy equa¬ 
tion of state w, the exact solution of Eq. Q for the grow¬ 
ing mode, neglecting the dark energy perturbations, is 
given by [BHH] 


5[a) = a 2 F 1 


1 1 
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for H{af/Hl = D.mo 
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where 2 Fiia,b;c] z) is a hypergeometric function, see 
Ref. [5] for more details. In more general cases, for in¬ 
stance admitting that the dark energy equation of state 
parameter is a function of time, it is impossible to find a 
closed form analytical solution for Eq. 0. 

However, when we take into account the effect of the 
dark energy perturbations, even though the effect is 
small, the analytical solution 0 is no longer valid. The 
reason for this is that, as shown in Ref. [S], the dark en¬ 
ergy perturbations can be included effectively as a Geff 
in Eq. 0: given by: 

1 — n 

Geff(«, k)/G^ = 1 + ^^(1 + . 

(3) 

Depending on the values of the equation of state parame¬ 
ter w and sound speed cl the density contrast of the dark 
energy can be of the order of A^e ~ 0.06Am and of the 
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order of Ade ~ O.OlAm, for w = —0.8 and w = —0.95 
respectively (of course if w = — 1 we are dealing with a 
cosmological constant and hence, by definition it has no 
perturbations). 

As mentioned before, the dark energy fluid considered 
in this paper should be thought as an effective dark en¬ 
ergy component, so the values of the equation of state w 
and sound speed are effective parameters, hence not 
always they have a physical meaning but rather they are 
a representation of the particular model taken into con¬ 
sideration. In practice, considering for instance a value 
of = 0 does not necessarily means that we are deal¬ 
ing with a dark energy fluid that behaves like dark mat¬ 
ter but it might be due to a particular modified gravity 
model expressed in terms of fluid parameters which has 
a value of cl equal to zero (or very small), see [TD] . It is 
also worth mentioning that the Eq. ^ has been evalu¬ 
ated under the condition of constant in time sound speed, 
however the sound speed can depend on the scale. 

The paper is organized as followed: in Section [l^ we 
find exact analytical expressions for the evolution of the 
matter density contrast sourced by the dark energy per¬ 
turbations and in Section |IV| we compare them to nu¬ 


merical solutions and test their accuracy; in Section III 


we find the new growth index 7 and we compare it with 
other parameterizations, while in Section |VI| we present 
our conclusions. 


II. THE SOLUTION FOR THE MATTER 
DENSITY CONTRAST 


In what follows we will find the analytic solution to the 
matter density contrast when dark energy perturbation 
are properly taken into account. To solve the second 
order differential equation for Sm we will consider two 
different limits: first when the sound speed is effectively 
zero {cl = 0 ), and second when the sound speed is 0 < 
cl < 1 . 

We start by rewriting Eq. Q in terms of G{a) de¬ 
fined as 6 {a) = aG{a), the dimensionless Hubble pa¬ 
rameter E{a) = H{z)/Hq and the parameter Q{a,k) = 
Ges{a,k)/GN, so Eq. Q becomes 


a‘^G"{a) + 


5 -I- a 


3 -t- a 


E'ja) 

E{a) 
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aG'{a) 


E'{a) _ 

E{a) 2 a^E{aY 


Q{g 


G(a)=0. (4) 


We need to express the function Q(a, k) which accounts 
for the dark energy perturbations^. In this paper we use 
the expressions found in m- In the latter the authors 
solved analytically the full system of differential equa¬ 
tions for dark matter and dark energy. In order to find 


^ we would like to remind the reader the G^f f and Q are the same, 
hence we use them interchangeably 


analytic solutions to matter and dark energy density con¬ 
trasts, some assumptions had to be made, in particular 
that the equation of state parameter w and the sound 
speed of the dark energy component have to be constant 
in time or at least slowly varying (however the sound 
speed can depend on the scale k). Once the analytic so¬ 
lutions for matter and dark energy were found then it was 
possible to express the quantity Q{a, k) which is defined 
as the ratio 


Q{a, k) = 1 + 


PdeAde 

PmAjyi 


( 5 ) 


which is a phenomenological function that takes into ac¬ 
count the relative growth of the dark energy perturba¬ 
tions. This function enters directly into the gravitational 
potential k'^cj) = —4:TTGa^Q{a, k)pmAm, that is the reason 
why Q{a, k) enters only in third term into Eq. Q. 

Another interesting issue is the initial conditions for 
matter and dark energy. In our case dark energy per¬ 
turbations are sourced by the dark matter perturbations 
via the gravitational potential; because we are in linear 
order perturbation theory then all fc-modes evolve inde¬ 
pendently and as a consequence each A:-mode depends 
linearly on a normalization factor which is constant in 
time however it depends on k, say 6o{k). The value of 
this constant is given by inflation at very early times. In 
this scenario, setting the initial conditions means to set 
the constant Jg. However, as stated, the factor Sg enters 
as multiplicative factor to both dark matter and dark 
energy, hence the quantity Q{a,k), which is the main 
concern of this paper, has no direct dependence on 
A detailed analysis on the initial conditions for the dark 
energy perturbations and the decaying modes of the so¬ 
lutions can be found in [TT]. The general solution, see 
Ref. [TT], is: 


Q{a,k) = 1-^ 1 ^mo (x + u,)- 

i Sr] 


where cl is the sound speed of dark energy. 


-, — 3w 
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A. The case = 0 


When the dark energy sound speed is equal to zero, 
then the Q{a, k) parameter which gives us the amount of 
the dark energy perturbations can be written as 

Q = 1 + Qoa-^^ ( 6 ) 

where 


_ 1 - Hmp 1+w 
Hmo 1 - 3w ■ 


( 7 ) 


Before inserting the above equation into Eq. Q, we can 
make the following change of variables to simplify the 
problem even more. Let us consider the new variable 


1 - n 


u = 
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then the Eq. Q becomes 

u^G + 

_1 1 - w - go 

1 + u 


^ 5 1 u 
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6w 21+ u 


uG = 0 




(9) 


where Qo = (1 + ~ Src). Eq. (H) can be solved 

analytically and the solution is 


G[a) — 2 F 1 

1 - a 


4 12w ’4 12w ’ 6w’ 
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( 10 ) 


where 




B. The case 7 ^ 0 

If the dark energy sound speed is not zero, then the 
dark energy will have a sound horizon below which per¬ 
turbations cannot grow. The modification to the Newto¬ 
nian potential will be, see m- 


Q{a) - 1 = 


O ( + ^) 2k^cla 


~ (1 - 11^0)0 ( 1 + W') 


tt2 

-^0 ^-l-3io 


^2^2 


, ( 12 ) 


where we have used the fact that we work in the sub¬ 
horizon approximation [k ^ aH). As it can be seen, 
the Q[a) term is suppressed by a term a~^ the lack of 
growth of the dark energy perturbations when they en¬ 
ter the sound horizon. However, inserting Eq. (12) into 
Eq. 0 and making the change of variable a —> m, we 
cannot find an analytic solution to the matter density 
contrast because of the extra dependence a~^ in the ex¬ 
pression of Q{a,k). In order to solve Eq. @ we make 
the approximation: 


Q{a) = 1-1- Qoa 

\3 


-l-3iu 


— I -|- Qqci 

H, 


— 3w 


(13) 


where Qo = (1 — Umo)! (1 + w) then the differen¬ 
tial equation that we need to solve will look exactly as 
Eq. (121 where now we have a new qq for modes below the 


r 2 

sound horizon: qo = ^ {1 + w) ° 2 fc”° ■ In this case we are 
overestimating the amount of dark energy perturbations 
below the sound horizon as we are taking out the term 
a~^ which lowers the clustering of the perturbations. 


Then the solution to the Eq. (|I2|) is 

Gia) 


2 El 
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where 




(15) 


C. The growth-rate 

In the previous section we have found the solution to 
the matter density contrast in terms of the hypergeomet¬ 
ric function: 


5{a) = 0 2^1 

I - 

where 

B = 


4 I2u> 4 1210 ’ 6w 


I 

flm(a) 


^V'(l-3u;)2 +24(55 (16) 

12 w 


and the parameter 6B can accommodate all different 
cases if defined as 


6B = 0 


SB = 


(1 -I- w) 
1 — 3?c 


(No DE perts), (17) 
(for cl = 0), (18) 




Then, the growth rate /(a) is 


(for cl > 0). (19) 


/(„) = 


+ a- 
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where the coefficients are 


2 El 




7 = -+a + (3. 

It is interesting to notice that the hypergeometric func¬ 
tion can be simplified making use of the relation be¬ 
tween hypergeometric functions and Legendre polyno¬ 
mial PI3{x). Using Eqs. (15.4.12) and (15.4.21) and 
Eqs. (8.5.3) and (8.5.5) of Ref. [S] we find 


/(a) = 1-1- 3wa I 1 - \/Hm(a)- 


jj5/6w 

2B-Si 

l/\/^2m(a) 

p5/6i(j 
25-i 
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If dark energy perturbations are switched to zero we have 


III. THE GROWTH INDEX 


/(a) = A/rJm(a)- 


p5/6w 

^ll6w 


p5/6w 
-1 /6w 

l/\/Hi„(a) 


( 21 ) 


in agreement with the expression found in Ref. 


D. Joint solution 

Alternatively, we can also find a joint solution to the 
equations in order to have one single solution. The dif¬ 
ference in the two solutions comes in the coefficients 
of the hypergeometric functions which then account for 
the dark energy perturbations, i.e. the coefficients B in 
Eqs. 0 and ( |15[ ). Hence we can think to join these two 
coefficients directly. We found that 


Sc2>o = Y^y(l-3u>)2-|-36Hi„o(l-l-'ic) 
y/il - + 6Bi 


Hi 

k'^cl 


B. 


c?=0 — 


12 u> 


( 22 ) 


(23) 


for scales above and below the sound horizon, respec¬ 
tively. Joining only the parts that account for dark en¬ 
ergy perturbations we find 


B 
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6 B 16 B 2 


dB,+SB2 


= Tb-J(l-M'+24- 

12w\l 1 _ + 2 


fe 2 


(24) 
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We can also add an a dependence to the coefficient. The 
reason is that the second term in the square root is ex¬ 
actly Q used, hence the dark energy perturbations, and 
we have 


Bjoint — ^ 2 ^ \/ (1 iwY + 


SB 1 SB 2 
5Bi + SB 2 


= ^,/(l-3w^)" + 24-- 

12«^V l-3u; + i^^c2 


.(25) 


However, in this paper we consider Eq. (24) for two main 


reasons: the presence of the scale factor does not change 
the solution for fas (the difference between the two is of 
the order of 10 “®, and second it is not really clear how 
to deal, in the hypergeometric function, with coefficients 
that depend on the variable. 

The expression found will be very useful for instance 
to speed up the code for forecasts, analytic estimates of 
7 or fitting the data. The growth index 7 will be shown 
in the next section to be the same, except that we have 
to use now the new 6B = 24- .2 -. 


Here we present the corrections to the growth index 
7 ( 0 ) due to the dark energy perturbations. In Ref. [12] it 
was shown that the growth rate /(a) = H-, in the case 
of no DE perturbations, can be approximated as 


/(a) = H„(a)^(“) 


^moia) = 


7(a) = 






-,-3 


H{aYim 

In/(a) 


3(1 — w) 
InHm(a) 5 — 6w 


(26) 

(27) 

(28) 


When we want to include the DE perturbations, one way 
to do it is via the semianalytic approach of Ref. m, 
where it was shown that the growth index depends on 


7 = 


3(1 - n; - A(Geff)) 


5 — 6w 

1 S 2 ni(Q') 


(29) 

(30) 


For other approaches that explore the effects of the dark 
energy perturbations or modified gravity in general on 
the growth index, see Ref. [TM2(1] . However, here we 
will follow a more straight-forward approach by using 
the analytic expressions found in the previous sections. 
In what follows we will use the shorthand of H to mean 
(a) in order to simplify the notation. Then using the 
ansatz /(H) = we find that the growth index up 

to first order can be written as 


7 


_ ln(/(H)) 
ln(H) 

3(JR-bu>- 1) 


6w — 5 

3(H - 1)((JR + w- l){9SB{4w - 3) - 3w + 2)) 
2((5-6w)2(12w-5)) 


-b ■ 


(31) 


The reason we perform the series expansion of 7 in terms 
of H is that it allows us to extract the zero-th order part 
of 7 , which in the ACDM model is 6/11, but also the 
first order correction. As a result we can compare with 
the other expressions commonly found in the literature, 
eg Eq. (28), and compare their accuracy. 


It is instructive to split the contributions to the growth 
index to two parts: 


7 = 7m +1DE, 


(32) 


where 


1m — 


3iw - 1) 3(3n; - 2)(w - 1)(H - 1) 

6w — 5 
3SB 


2(5-6w)2(12w-5) 

3SB{6w{6w - 11) -b 29) 


, (33) 


'2((5-6r(;)2(12w-5)) 


27SB'^{Aw - 3) 
2((5-6w)2(12w-5)) 
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(34) 












































5 


where the first term 7 ^ corresponds to the well known re¬ 
sult from Ref. [12] and the second term 'jde corresponds 
to the extra contribution of the DE perturbations. 

Alternatively, we can follow the procedure of Ref. [T^ 
to evaluate the growth rate. We can do this by changing 
variables from a to and replacing the growth factor 
5{a) with the growth rate /(fl) = Doing that we 

arrive at a differential equation for /(fl) in terms of fl: 

3wii - n)nf'(n) + f{n) Q - ^w{i - + 

= 0 . (35) 

Using the ansatz /(U) = we arrive to the same 

result as in Eq. (j^. 


A. The case cj = 0 



Figure 1: The percent difference of the effect of the DE per¬ 
turbations of 7 as part of the total and as a function of the 
equation of state w for = 0. In practice we plot 100 * 
with the latter given by the zero order part of Eq. (40l. 


In the case when the sound speed is zero, we can use 
Eq. (32) with SB = 73 ^ to find 


7(D) 


3w{3w — 5) 

{3w — l)( 6 r(; — 5) 

15r(;(3w - 5) - 5) 

2(12r(; - 5) (18^2 - 21w-f 5)^ 

(36) 


Also, we separate the effects of the matter and dark 
energy perturbations on 7 as 


7 = 7m + iDE, (37) 

where 

_ 3(u;-l) 3(3u;-2)(u;-l)(U-l) 

6w-5 2(5-6w)2(12ri;-5) ’ ^ ^ 

_ 3{w + l) 

18r(;2 _ 21u; -b 5 

in fx 3(w-b l)(9w(3r(;(2w — 5 )-b 8) — 1 ) 
(12u;-5)(18u;2_21i(;-b5)^ 

(39) 

We can also compute the ratio of the two quantities 
iDE/lm in order to estimate the effect of neglecting the 
DE perturbations on the measurement of 7 in up-coming 
surveys. We find that the ratio can be written as 

iDE _ W + l ^ 

7m 3r(;2 — 4 ^; -b 1 

9w{w + l){3w — 5)(4'u; — 3)(fl — 1) 

2(1 - 3w)2(u> - l)(6w - 5)(12w - 5) ’ ■ 

(40) 

In Fig. [2 we show the percent difference of the ef¬ 
fect of the DE perturbations as part of the total and as 
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Figure 2: A comparison of all the different expressions for 7 
for (c2,w,f2mo) = (0,—0.8, 0.3). The solid black line corre¬ 
sponds to the expression for the analytic solution of Eq. (36 1 , 


the dashed black line to the numerical solution from CAMB, 
the dotted line to the wCDM model with no DE perturba¬ 
tions, the dot-dashed black line to the semianalytic expression 
of Eq. (291, while the red and dashed-red lines correspond to 


the zero-order terms for expansions for 7 without and with 
DE perturbations respectively. 


a function of the equation of state w. In practice we 
plot 100 * with the latter given by Eq. (40). For 
w = — 1 Eq. (36) gives the expected result 7 = 6/11, 
but for w ^ —1 there are substantial differences, eg for 
w = —0.8 we have 7 = 0.533 which is a 3.4% differ¬ 
ence from the result 7 = 0.551 found when neglecting 
DE perturbations. These differences can be important 
as cosmology has moved into a high precision era and 
percent accuracies will be sought after in the cosmologi¬ 
cal parameters by the upcoming surveys like DES, LSST 
and Euclid, where the relative errors on the growth in¬ 
dex 7 are of the order of 0.5 — 3% depending if combined 
probes are used, see 0 , m, m, m- 
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In Fig. [^we show a comparison of all the different ex¬ 
pressions for 7 for fimo) = (0, —0.8,0.3). The solid 

black line corresponds to the expression for the analytic 
solution of Eq. (36), the dashed black line to the numer¬ 
ical solution from CAMB, the dotted line to the wCDM 
model with no DE perturbations, the dot-dashed black 
line to the semianalytic expression of Eq. (29), while the 
red and dashed-red lines correspond to the zero-order 
terms for expansions for 7 without and with DE pertur¬ 
bations respectively. 

Clearly, our analytic solution in this case is in excellent 
agreement with the numerical one and by far superior 
to the semianalytic one. Specifically, the expression of 
Eq. (29), which is the dot-dashed curve in Fig. differs 
by more than ~ 6 % at z ^ 0. The reason for this is 
that it is a semianalytic approximation derived at high 2 ; 
and as a result at low redshifts it is significantly different 
from both the exact and numerical solutions. Therefore, 
we find that it is not suitable for use in Fisher Matrix 
analysis or data fitting when the DE perturbations are 
taken into account. 

Einally, we should mention that we have checked that 
using the growth index 7 of Eq. © to calculate the 
growth-rate / and fa% by using and integrating Eq. (26) 
respectively, is in agreement to better than 0 . 1 % com¬ 
pared to using the full analytical solution of Eq. (10). 


B. The case 7 ^ 0 


When the sound speed is different from zero, we can 
again use Eq. (32) with 6B = ||k( 1 -l-w), where we have 


set K = 




k'^c 


to find 


_ , 3(3k- b (3 k- b 2)r(; - 2) 3(0-1) 

’ ~ 2{Qw - 5) 8(5-6i(;)2(12u;-5) ' 

(27k + w — 3) — Ore -b 4) • 

(3K-b(3K-b2)w-2)-b--- . (41) 


As mentioned earlier in this section and in order to sim¬ 
plify the notation, by writing O we actually mean (®) 
and Om,o = Oi„o(a = 1 ). 

Also, we can separate the effects of the matter and 
dark energy perturbations on 7 as 

7 = 7m + IDE, (42) 

where 

3(w - 1) 3(3w-2)(w- 1)(0- 1) 

6w-5 2(5-6w)2(12w-5) ’ ^ ^ 

9li(w + 1) 9(0 — l)K(?i' -b 1) 

^ 2(6w- 5) ~ 8(5-6r(;)2(12r(;-5) ' 

(27k (4r(;2 -|- _ 3) q- I2w{6w — 11) -b 58) . 

(44) 


We can also compute the ratio of the two quantities 
iDE/lm in order to estimate the effect of neglecting the 


DE perturbations on the measurement of 7 in up-coming 
surveys. We find that the ratio can be written as 


7£>g _ 3 k(w -b 1) 

7m 2(w - 1) 

(4u> — 3) (3 k 


- 27(0 - 1)k(u;- b 1) 
(3 k -b 2)w — 2) 


8 (w- l)(6w-5)(12w-5) 


(45) 


IV. COMPARISON WITH THE NUMERICAL 
SOLUTION 


Here we compare our analytical solution of Eq. (14) 
with the parameter B given by Eqs. (11) and Eqs. (15), 
with the numerical one for various choices of parameters 
and for both cases. Since the data are given in terms of 
the parameter fas{a) = f{a)a{a) = |^ad'(a), we will 

prefer to compare that and the growth rate /(a) = q/p® ■ 

In Fig. (left) we show a comparison of the analyt¬ 
ical solution for f(j8{z) (dashed line) vs the numerical 
solution of the full system of differential equations for 
the DE perturbations by CAMB (solid black line), the 
case with no DE perturbations (dotted line) and the 
f2m(a)''' parametrization with 7 given by Eq. (29) (dot- 
dashed line) for = 0 , fc = 20077o and (w, Umo: cr 8 ,o) = 
(—0.8,0.3, 0.8). The gray points correspond to mock /cr 8 
data based on the specifications of an LSST-like survey 
m, as it was done in Ref. m- As can be seen, the differ¬ 
ence between the cases where we include (black solid or 
dashed lines) or neglect DE perturbations (dotted line) 
can be significant compared to the small error bars of the 
expected data from Euclid or LSST and might bias the 
results. 

In Eig. ^ (right) we show of the analytical solution 
(black dashed line) vs the numerical solution (solid black 
line) of the ODE of Eq. ([^ for = 1 and k = IOHq. 
In this case, we also test the effect of the different scale 
k but also the non-zero sound speed and we see that 
the difference can be smaller. Our analytical solutions 
either exact or approximate were found to be in excellent 
agreement with the numerical ones. In Fig. we show 
the percent difference between the analytical solution and 
the numerical one from CAMB, for various values of the 
sound speed and the scale k and as it can be seen, in 
all cases we find agreement better than 0.5%. 


V. FORECASTS 


In this section we test the new expression for the 
growth index 7 given by Eq. (31) in order to constrain 


the cosmological parameters with the upcoming surveys 
and we compare the results with the constraints ob¬ 
tained using the growth index given by Eq. (29). To 


constrain the cosmological parameters such as the equa¬ 
tion of state parameter w and the sound speed Cg, we 
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Figure 3: Left: Comparison of the analytical solution for fa8{z) (dashed line) vs the numerical solution of the full system 
of differential equations for the DE perturbations by CA MB (solid black line), the case with no DE perturbations (dotted 
line) and the flm(a)^ parametrization with 7 given by Eq. (291 (dot-dashed line) for = 0, fc = 200770 and (w, ffmo, rrg.o) = 
(—0.8,0.3,0.8). The gray points correspond to mock fa8 data based on the specifications of an LSST-like survey m , as it was 
done in Ref. [13]. Right: Comparison of the analytical solution (black dashed line) vs the numerical solution (solid black line) 


of the ODE of Eq. (^Ik for c^ = 1 and k = IOHq. 
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Figure 4: The percent difference between the analytical solu¬ 
tion and the numerical one from CAMB, for various values of 
the sound speed and the scale k. 


adopt the Fisher matrix technique having in mind a set 
up similar to LSST experiment [22|. The LSST is ca¬ 
pable of exploring the universe in a range in redshifts of 
2 ; € [0,2] covering an area of about 20000 deg^. The er¬ 
ror on redshift is Uz = 0 . 02(1 -|- z) which corresponds to 
the goal of the LSST; the galaxies are distributed with 
n{z) oc exp [—(z/zi)^] where a = 2, Zi = 0.5 and 

13 = 1 . 

In order to avoid non-linearity problems (both in the 
spectrum and in the bias), we evaluate the Fisher ma¬ 
trix up to a limiting kmax[z) at each z: we choose 
values from O.llh/Mpc for low-z bins to 0.3/i/Mpc 
for the highest z-bins, see Ref. [231125] for more de¬ 
tails on the Fisher matrix calculations. For the Fisher 
matrix we consider the following parameters: 9 = 
[uJni,UJhjT,ns,^mfi,W,cl,Pshot], where UJm = 


ujh = and Pshot the shot noise. The reference cos¬ 

mology used in the analysis is the best-fit given by Planck 
m with the exception of w = — 0.8 and the sound speed 
Cj for which we used first a value = 0 and then = 1 . 

For a sound speed of = 0 we find: {ctujjO’cs} = 
{0.0089,1.67 X 10"®} and {cr^,crc2} = (0.0104,3.69 x 
10“®}, for that of Eq. (29) and our expression respec¬ 
tively. For a sound speed = 1 we find: {aw,<Jc'^} = 
{0.0052,76.63} and = {0.0058,94.00}, for the 

semianalytic expression and the analytic expression re¬ 
spectively. 

Using the expression for the growth index given by 
Eq. (31), the errors on the parameters are increased of 
about 50% for the sound speed and of about 15% for the 
equation of state parameter, for the = 0 case. The 
reason can be found looking at the Fig. The semi¬ 
analytic expression for 7 (a) especially at low redshifts is 
smaller than the full numerical result: the growth rate 
/(a) (and also the growth of matter G{a)) are propor¬ 
tional to and so do the derivatives with respect 

to the cosmological parameters. The matter density pa¬ 
rameter Um(a) is always lower the unity (it is equal to 
1 only at very high redshifts during matter domination 
era), consequently if the growth index 7 is decreased the 
overall effect on the growth rate is increased and this is 
reflected in Fig. Hence, using the semianalytic expres¬ 
sion overestimates (i.e. reduce) the errors. 


VI. CONCLUSIONS 

In Ref. |5] it was shown that the effect of the dark 
energy perturbations on the matter density perturbations 
is to induce a Geff(a, fc), as shown in Eqs. 0 and 0. 
In this work we have found the analytical solution to 


















Eq. Q for constant w when = 0 and we have explicitly 
shown that not only is the solution different from the one 
when the DE perturbations are neglected, see Eq. ([^, 
but the difference is actually large enough to affect the 
results of future surveys like Euclid or LSST, see Fig. 
(right). In the case of 0 < cf < 1 we presented analytic 
approximations that are in excellent agreement with the 
numerical solution to Eq. Q. 

Furthermore, we also found analytical expressions for 
the growth index 7 used commonly to parameterize the 
growth as /(a) = . This parameter is given by 

7 ~ ^ model with constant w in GR when 

the DE perturbations are ignored and is 6/11 for ACDM 
(w = —1), but when the DE perturbations are properly 
taken into account it is given by 7 ~ ( 3 ^-i^)™ 6 ^^ 25 ) when 
Cg = 0 , which however again becomes 6/11 for w = —1 
since ACDM does not have DE perturbations. 

In the case when the DE sound speed is non-zero 
we found that the growth index is given by 7 ~ 
3 ( 3 K+( 3 K+ 2 )m- 2 ) ^ ^ and has an explicit 

scale dependence. Our expression is different but at the 
same time more accurate than other expressions that 
have appeared in the literature. We found that it works 
to better than 0.3% accuracy for k > IOHq or equiv¬ 
alently k/h > 0.0033Mpc“^, thus making it extremely 


useful for use in forecasts for future surveys. We have also 
compared our new expressions to that one of Eq. (291, 
via a Fisher Matrix approach, and we found that the lat¬ 
ter significantly overestimates the errors compared to the 
more accurate growth index of Eq. (31) we presented for 
the first time in this paper. 

In conclusion, we have shown that using less accu¬ 
rate expressions for the growth index, such as the one 
of Eq. (29), or more importantly neglecting completely 
the DE perturbations can lead to misleading estimations 
of the growth index 7 to the percent level, something 
which is highly relevant as cosmology today has moved 
into a high precision era and percent accuracies will be 
sought after in the next generation surveys. 
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